Bose-Einstein condensation of trapped atoms with dipole interactions 
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The path integral Monte Carlo method is used to simulate dilute trapped Bose gases and to inves- 
tigate the equilibrium properties at finite temperatures. The quantum particles have a long-range 
dipole-dipole interaction and a short-range s-wave interaction. Using an anisotropic pseudopoten- 
tial for the long-range dipolar interaction and a hard-sphere potential for the short-range s-wave 
interaction, we calculate the energetics and structural properties as a function of temperature and 
the number of particles. Also, in order to determine the effects of dipole-dipole forces and the 
influence of the trapping field on the dipolar condensate, we use two cylindrically symmetric har- 
monic confinements (a cigar-shaped trap and a disk-shaped trap). We find that the net effect of 
dipole-dipole interactions is governed by the trapping geometry. For a cigar-shaped trap, the net 
contribution of dipolar interactions is attractive and the shrinking of the density profiles is observed. 
For a disk-shaped trap, the net effect of long-range dipolar forces is repulsive and the density profiles 
expand. 

PACS numbers: 03.75.Hh, 03.75.Nt, 05.30.Jp, 02.70.Ss 



CO 



I 

C 

O 

o 



> 

00 

^t- 

ON 

o 
in 
o 



-a 
c 

o 
o 



X 



I. INTRODUCTION 

The experimental observation of Bose-Einstein con- 
densation (BEC) in dilute trapped and supercooled 
atomic vapors |l| has stimulated enormous experimen- 
tal and theoretical interest in the physics of the weakly 
interacting Bose gas (see 2] for a review). Especially in- 
triguing is the role of dimensionality and two-body inter- 
actions. Lower dimensional atomic gases have been real- 
ized with strong quantum confinement in one or more di- 
rections. These systems exhibit very peculiar properties, 
such as the exact mapping between interacting bosons 
and noninteracting fermions (Tonks-Girardeau gas) 
because the role of quantum fluctuations and strong cor- 
relations is enhanced. Although these gases are very 
dilute, BEC of trapped atomic gases is strongly influ- 
enced by two-body interactions. Within a mean-field 
approximation, the strength of the interatomic interac- 
tions can be typically characterized by a single parame- 
ter, the s-wave scattering length a s , at very low temper- 
atures. This van der Waals type interaction is isotropic 
and short-range in comparison to the average interpar- 
ticle separation. Accordingly, Bose-Einstein condensates 
created thus far have been very accurately described by 
the Gross-Pitaevskii theory (see Q for a review) , a mean- 
field approach. Furthermore, it has been possible to tune 
the s-wave scattering length a s to essentially any value, 
positive or negative, by utilizing a magnetic atom-atom 
Feshbach resonance Q] . 

Studies of the influence of the interparticle interac- 
tions on the properties of BEC have focused mainly on 
the short-range interactions Q. However, understand- 
ing new kinds of systems will require the investigation of 
the effects of additional interactions, since Bose-Einstein 
condensation of molecules has now been observed using 
new techniques for creating ultracold molecules In 
addition, Ling et al. Q proposed an efficient method to 
convert atomic condensates into stable molecules using a 



generalized Raman adiabatic passage scheme. Very re- 
cently, the generation of a BEC in a gas of 52 Cr atoms 
was reported Q. As molecules can potentially have 
large dipole moments and a chromium atom has a very 
high magnetic moment of 6/ib (Bohr magneton) in its 
ground state, BEC with dipole-dipole interactions in ul- 
tracold gases has attracted considerable theoretical at- 
tention over the past few years H, 0, 0, 0, 0, 0, 
Q E E El Compared to the short -range in- 

teraction, the dipole-dipole interactions are long-range, 
anisotropic, and partially attractive. As a result, nonlin- 
ear effects due to long-range dipole-dipole interactions 
lead to new interesting properties; different quantum 
phases such as superfluid, supersolid, Mott insulator, 
checkerboard, and collapse phases E3> a self-bound Bose 
condensate in the field of a traveling wave Jlal . a ro- 
ton minimum in the excitation spectrum [Tj, Il5| . and 
promising candidates for the implementation of fast and 
robust quantum-computing schemes E3|- First of all, Yi 
and You suggested an effective pseudopotential for 
anisotropic dipolar interactions in the Born approxima- 
tion and recently Derevianko |To) derived a more rigorous 
velocity-dependent pseudopotential for dipole-dipole in- 
teractions. Santos et al. 11] showed a strong dependence 
of the stability of trapped dipolar BEC on the trapping 
geometry. O'Dell et al. |l7j derived exact results in the 
Thomas-Fermi regime for the statics and dynamics of 
a trapped dipolar BEC. However, in all the theoretical 
studies so far the T=0 K situation has been investigated. 
In this paper we address the properties of trapped con- 
densates with dipolar interactions at finite temperatures. 
Several sources for dipolar BEC are atoms polarized by 
an electric field or atoms with large magnetic moments 
in addition to polar molecules. The magnitude and sign 
of dipole-dipole interactions can be changed by rapidly 
rotating an external field |l8j . Thus, by varying the 
strength of the dipolar interactions relative to the s-wave 
scattering (tuning the dipolar interaction or changing the 
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s-wave scattering length a s ) it is possible to study the ef- 
fects of the dipolar interaction on a condensate in detail. 

In this paper we attempt to understand Bose-Einstein 
condensation of trapped dipolar gases, where the inter- 
particle interaction consists of a long-range dipole-dipole 
interaction in addition to the usual short-range s-wave 
interaction (a purely repulsive hard-sphere potential). In 
particular, we demonstrate the influence of dipole-dipole 
forces on the properties of ultracold gases by compar- 
ing the energetics and structural properties of harmon- 
ically trapped atomic condensates with dipolar interac- 
tions with those for which only short-range interactions 
are present. In addition, we study the dependence of the 
sign and the value of the dipole-dipole interaction energy 
on the trapping geometry by using two different types of 
cylindrical traps ( a cigar-shaped trap and a disk-shaped 
trap). For these purposes, we use a finite-temperature 
path-integral Monte Carlo (PIMC) method |Ts| . 

Our paper is structured as follow. SectionlTTIintroduces 
the Hamiltonian of trapped dipolar hard-sphere systems 
used in this studies and an effective pseudopotential to 
describe anisotropic dipolar interactions in the Born ap- 
proximation. The path integral Monte Carlo technique 
briefly is also reviewed in Sec. [HI Section lrTfl presents our 
simulation results. Finally, we conclude in Sec. IIVI 



II. DESCRIPTION OF THE SYSTEM AND 
SIMULATION METHOD 

We wish to study the problem of a quantum iV-particle 
system to consider a harmonically trapped Bose-Einstein 
condensation in three dimensions. The quantum parti- 
cles have both a long-range dipole-dipole interaction and 
the usual short-range hard-sphere interaction. We as- 
sume that all atomic dipole moments are equal and ori- 
ented along the z-axis (trap axis). Within the Born ap- 
proximation, the dipole-dipole interaction potential be- 
tween two dipoles separated by is given by 







(r > a s ). 



(3) 



/i 2 (l - 3 cos 2 9ij) 



(1) 



where p characterizes the dipole moment, r.y = rj — rj 
is the vector between the dipoles i and j, and 9ij the an- 
gle between the vector ry and the direction of the dipoles 
(z). To date, this pseudopotential has worked reasonably 
well for the long-range dipole-dipole interaction. The cor- 
responding Hamiltonian is then 
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H = H + J2v d (n 3 ) + J2 v ( r v) ( 2 ) 

i>j i>j 

h 2 N l N 



i=l i=l 

where Hq is the Hamiltonian for trapped ideal boson 
gases and v(r) is the hard-sphere potential defined by 

v(r) = +oo (r < a s ) 



In order to investigate the dependence of the sign and 
the value of the dipole-dipole interaction energy on the 
trapping geometry, we consider Bose gases under a cylin- 
drically harmonic confinement (oj x — oj v — Xlo z and u> z ), 
where A is the aspect ratio. 

Consider a system of N dipolar hard spheres. The 
statistical mechanics of a quantum system at an inverse 
temperature /3, where = (ksT)^ 1 and fc^ denotes the 
Boltzmann constant, is governed by the density matrix. 
The Bose-symmetrized density matrix is given by 



(4) 



where 1Z and TV are two configurations of N dipolar 
hard spheres, which denote 3N dimensional vectors, 1Z 
= (R-1,- ■ -,Rat). P denotes a permutation of particle 
labels among hard spheres and P1Z is one such permu- 
tation. The inclusion of particle permutations is crucial 
and superfluidity is an immediate consequence of Bose 
symmetry. However, evaluating the density matrix for 
interacting systems at very low temperatures is compli- 
cated by the fact that the kinetic and potential terms in 
the exponent of the density matrix do not commute. In 
order to avoid this problem, we insert M—l intermediate 
configurations into Eq. Jj}: 

p(K, PW; 13) = ...J dKidK 2 ■ ■ ■ dU M -i 

x p (H,Ki;t) ■ ■ ■ p(TZ M -i,Pn';T),(5) 

where r = [3/M is the imaginary time step and M is 
referred to as the number of time slices. Eq. [S] is the 
path-integral formulation of the density matrix, which 
contains the discrete sum over permutations and multi- 
ple integrations over density matrices at a higher tem- 
perature (fcer) -1 . Eq. can be evaluated in the path 
integral Monte Carlo technique by a stochastic sampling 
of the discrete paths {7?., 7^i, 7^2, • ■ ■■7£m-i, P1Z'} using 
multilevel Monte Carlo sampling |19ll20l |. which accounts 
for permutations. In order to use Monte Carlo sampling, 
we must first provide a pair-product form of the exact 
two-body density matrices. For the short-range hard- 
sphere potential, we used the high-temperature approx- 
imation for the hard-sphere propagator derived by Cao 
and Berne pll . Unfortunately, it is impractical to de- 
velop an analytic pair-product form for the anisotropic 
dipole-dipole interaction potential. Thus the primitive 
approximation, which separates the density matrix into 
kinetic and potential energy terms, was used. In order 
to choose the value of the number of time slices M, we 
performed consistency checks by varying M to see that 
the results had converged. We used up to 160 time slices 
(M = 160) depending on the temperature and the num- 
ber of particles. We employed the canonical ensemble, 
i.e. in each simulation we fixed the temperature T, the 
number of particles N, the diameter of a hard sphere 
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a s , and the strength of the dipolar interactions \x. The 
thermal expectation value of an operator O is given by 
(O) = Tr(psO)/Tr(ps)- Expectation values are cal- 
culated at the same time the discrete paths are sam- 
pled. Typically 20,000 - 30,000 MC steps are required for 
equilibration and statistical averages are collected from 
150,000 - 180,000 MC steps after this. In each MC step, 
we attempted 60 - 1,250 trial moves at each time slice. 
Multiple runs were performed for each set of parameters 
so that statistical error bars could be determined. We 
performed these calculations partly on IBM Power4 and 
partly on Athlon processors at the University of Georgia. 
For N = 125 dipolar particles and r = 0.00446, the simu- 
lation (thermalization + computing averages) at a single 
temperature took about 2100 IBM Power4 CPU hours. 
The density profile n(R), which indicates the likelihood 
of finding particles at a given distance R from the center 
of the trap, is defined as 



n{R) 



M-l N 



M 



K= i=l 



S(R t , K - R) 
R 2 



(6) 



where Ri tK denotes the distance between the position of 
the ith particle and the center of the trap at time slice k 
and the symbol < • > denotes an ensemble average. The 
superfluid fraction along the axis of rotation n is given 
by0 
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p (3h 2 i° 

where A n — An denotes the projected area, 
M-l N 

2 



^ M-l N 
A = O H( Ri > K X R '^+l)' 



(7) 



(8) 



K=0 1 = 1 

and 1° the classical moment of inertia, 

M-l N 

M ' 



/ M-l N \ 

J ° = ( m E E TT( n x R i,«) • (n x R ijK+1 ) 



(9) 



3 0.34 




50 100 1 

Number of Particles N 



0.25 



0.05 




0.2 


N = 27 








N = 64 


profile 

p 


N = 125 


Density 

p 





1=0 i=l 



FIG. 1: The calculated density profiles in a cigar-shaped trap 
(a) n(p) and (b)n(z) at T = 0.4 T C 3D for three different num- 
bers of particles, N = 27, 64, and 125. The inset in (a) shows 
half widths. For comparison, we used two values fi = (open 
symbols) and n — 0.2 (filled symbols). A thick solid line in 
(a) is the density profile n(p) for an ideal Bose gas. In all 
figures, when statistical errors cannot be seen on the scale of 
the figure, the error bars are smaller than the symbol sizes. 



III. SIMULATION RESULTS AND DISCUSSION 

The path integral Monte Carlo (PIMC) method allows 
one to calculate accurate quantum mechanical expecta- 
tion values of many-body systems, the only input being 
the many-body potential. Here, we use PIMC to simulate 
a system composed of dipolar hard spheres for various 
values of the number of hard spheres N and the trap as- 
pect ratio A at two fixed dipolar moments \x as a function 
of temperature T and to obtain simulation results for a 
Bose-Einstein condensation with magnetic dipole-dipolc 
forces. This section discusses the energetics and struc- 
tural properties of trapped dipolar hard spheres obtained 
from our PIMC simulations. For simplicity, throughout 



this paper, we use a cylindrically symmetric harmonic 
trap with an aspect ratio A, A = lo x /uj z and lu x = uj y . 
The approximate critical temperature T^ D for an ideal 
gas of N atoms in a trap is given by [2( 

W 3 ) 2u ho [5 (3)] 2 / 3 W) ' ( j 

3D _ fiuho / N \ 1/3 

k B U(3)y ' 

where Tq D is the critical temperature for an ideal Bose 
gas in a trap in the thermodynamic limit, g(n) is the 
Riemann zeta function, and Co = (cu x + uj y + lu z )/3 and 
ujho = (ux^yi^z) 1 ^ 3 are the arithmetic average and the ge- 
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FIG. 2: In a cigar-shaped trap, the total energy per particle 
E/N as a function of the number of particles N at T — 0.4 
Tc° for fx = (open circles), p = 0.2 (filled circles), and an 
ideal Bose gas (open squares). 



FIG. 3: In a cigar-shaped trap, the dependence of the density 
profile n(p) on the temperature T for N = 64. 



ometric average of the trapping frequencies, respectively. 
A measure of the strength of the dipolar interaction rel- 
ative to the s-wave scattering in general is given by 



--dd 



p 2 m 
3h 2 a s 



(11) 



For chromium, the more common bosonic isotope, 52 Cr, 
has Edd ~ 0.089 and the less common bosonic isotope, 
50 Cr, has £dd ~ 0.36 [24[. In this study, in order that 
the effects of dipolar interactions can be made visible, 
we used m = m( 87 Rb), uj z {uj x ) = 2-rtx 77.78 Hz, uj x = 
10 lo z (lu z = 10 uj x ), the s-wave scattering length a s = 
10 a,Rb = 0.0433 a z (a x ), and the dipole moment p = 0.2 
(TuOzO^,) 1 ' 2 ((Hu) x a x ) 1 /' 2 ) for a cigar-shaped trap (a disk- 
shaped trap) as an example for numerical calculations, 
giving Edd ~ 0.31. Here, a^b — 100 a (Bohr radii) is the 
s-wave scattering length of 87 Rb, a z = ^h/muj z , and a x 
= y/U/muJx. For these parameters, the short-range in- 
teractions and the dipole-dipole ones have a comparable 
strength, and we can see the interplay between short- 
range and dipole-dipole interactions and the deviation 
between the case with and without dipole-dipole interac- 
tions. The results are scalable. Moreover, the magnitude 
and sign of the values of a s and fi can be easily tuned 
by Feshbach resonance Q and by an external field [l8[ . 
respectively. 



A. A cigar-shaped trap 

We first consider a 3D dipolar Bose gas in a cigar- 
shaped trap with A = tu x /uj z — 10 and lo x = uj y by in- 
creasing the confinement in the transverse direction. The 
length unit is a z = \JU/mu> z and energies are measured 
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FIG. 4: In a cigar-shaped trap, the superfluid fraction p s /p 
along the axis of rotation z as a function of the number of 
particles N at T = 0.4 T C 3D for fi = (open circles), fj, = 0.2 
(filled circles), and an ideal Bose gas (open squares). 



in units of hu z - Figure 1(a) shows the calculated den- 
sity profiles n(p) as a function of p, normalized such that 
J °° n{p)pdp = 1 and p = v / x 2 +y 2 , at T = 0.4 T C 3D for 
three different numbers of particles, N = 27 (circles), 64 
(diamonds), and 125 (squares). For comparison, we used 
two values p — and p = 0.2. Filled symbols indicate 
the density profiles n(p) for p = 0.2 and open symbols 
the density profiles n{p) for p = 0. In all figures, when 
statistical errors cannot be seen on the scale of the figure, 
the error bars are smaller than the symbol sizes. Qual- 
itatively, the density profiles for p — and for p = 0.2 
are similar in shape, as predicted. However, the density 
profiles n(fl) for p = are broader than those for p — 
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FIG. 5: The calculated density profiles in a disk-shaped trap 
(a) n(p) and (b) n{z) at T = 0.4 T C 3D for three different num- 
bers of particles, N = 27, 64, and 125. For comparison, wc 
used two values p = (open symbols) and p = 0.2 (filled 
symbols). A thick solid line in (b) is the density profile n(z) 
for an ideal Bose gas. The inset in (a) shows half widths. 



0.2. The difference with respect to the density profile 
n(p) without dipolar interactions is clearly seen as the 
number of particles increases. The density profiles n(z) 
along the axial direction as a function of z, normalized 
such that n(z)dz = 1, at T = 0.4 T* D for three dif- 
ferent numbers of particles are also shown in Fig. 1(b). 
In contrast to the density profile n(p), the difference in 
the density profiles n(z) between for p =0 (open symbols) 
and for p = 0.2 (filled symbols) is small as the number 
of particles grows. In addition, a thick solid line in Fig. 
1(a) shows the density profile n(p) for an ideal Bose gas 
for comparison. The density profile n(p) for the interact- 
ing gas approaches that for an ideal gas as the number of 
particles N decreases. Correspondingly the excitations 
in the longitudinal direction are largely frozen out and 
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FIG. 6: In a disk-shaped trap, the energy per particle E/N 
as a function of the number of particles N at T = 0.4 T% D 
for p — (open circles), p = 0.2 (filled circles), and an ideal 
Bose gas (open square). 



the system behaves quasi-one dimensionally. Our den- 
sity profile studies support the facts that i) as we add 
more particles into the gas, the density profiles expand; 
the half-width of the density profile increases with the 
number of particles N, as shown in the inset of Fig 1(a) 
and ii) the net effect of dipole-dipole interactions in the 
cigar-shaped trap is attractive, i.e. shrinking of the den- 
sity profiles is observed. To further characterize the be- 
havior of particles with dipolar interactions at T = 0.4 
T^ D , Fig. 2 shows the total energy per particle E/N of 
the system as a function of the number of particles N. 
The total energy per particle is strongly affected by the 
particle-particle interaction and the overall dipole-dipole 
interaction is attractive in the cigar-shaped trap, because 
the total energy per particle for p = 0.2 is smaller than 
that for p = 0. From Fig. 2 we can clearly see the de- 
pendence on the number of particles N, and numerical 
calculations have previously shown that the finite-size ef- 
fects are significant for rather small values of N [2^, [2j| . 
The difference in the total energy per particle E/N be- 
tween for an ideal Bose gas and for an interacting gas 
increases with increasing N, which is a consequence of 
the expansion of the density profiles with interaction, as 
shown in Fig. 1, although the ideal bosonic gas energy 
per particle E/N changes considerably less as the number 
of particles N increases. Accordingly, the dependence of 
the total energy per particle E/N on the number of parti- 
cles N reflects essentially the interaction energy alone, as 
the kinetic energy depends weakly on N . Figure 3 shows 
the dependence of the density profiles n(p) on the tem- 
perature T for N = 64. As the temperature increases, 
the effect of interactions is vanishingly small. At T = 
2.0 T^ D , the density profiles n(p) for p = and for p = 
0.2 are indistinguishable and the density profiles have the 
Gaussian profiles of a pure thermal distribution. Finally, 
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FIG. 7: In a disk-shaped trap, the superfluid fraction ps/p 
along the axis of rotation z as a function of the number of 
particles N at T = 0.4 T C 3D for p = (open circles), p = 0.2 
(filled circles), and an ideal Bose gas (open squares). 

we calculate the superfluid fraction p s j ' p along the axis of 
rotation z and in Fig. 4, we plot the superfluid fraction 
p s /p as a function of the number of particles N at T = 
0.4 T^ D . In addition, Fig. 4 shows the superfluid frac- 
tion p s /p for an ideal Bose gas. In general, the superfluid 
fraction is essentially unity at very low temperatures and 
decreases gradually with increasing temperature. For the 
ideal gas, the superfluid fractions show at most a very 
weak dependence on the number of particles considered. 
However, there are significant finite-size effects for the 
interacting gas. As the presence of repulsive interactions 
has the effect of lowering the value of the critical tem- 
perature, the superfluid fraction for the interacting gas 
is smaller than that for the ideal gas at the same temper- 
ature T/T C 3D , as shown in Fig. 4. The dependence of the 
superfluid fraction p s /p on the presence of interactions 
increases with the number of particles N. It is not easy, 
however, to see the difference between the superfluid frac- 
tions p s /p for p = 0.2 (filled circles) and those for p = 
(open circles) because of the size of the errorbars but 
the difference must clearly be small, if any. 



bosonic gas. Compared to the density profiles for A = 
10, they are similar in shape. As u) z is larger than u> x , 
the density profile n(z) for an interacting gas approaches 
that for an ideal gas (see Fig 5(b)), indicating the exci- 
tations in the axial direction are largely frozen out and 
the system behaves quasi-two dimensionally. In a disk- 
shaped trap, the density profiles expand along both the 
z-axis and the p-axis as we increase p., in agreement with 
those found by other authors. Correspondingly, the net 
contribution of long-range dipolar forces is repulsive, in 
contrast to the cigar-shaped trap. Simulations for other 
values of A show that BEC in a disk-shaped trap is sta- 
ble, even in very asymmetric traps, because the dipole- 
dipole interaction energy always remains positive. Figure 
6 shows the total energy per particle E/N as a function 
of the number of particles N both for interacting atoms 
with p, = (open circles) and p — 0.2 (filled circles) and 
for noninteracting particles (open squares). The total en- 
ergy per particle E/N for p — 0.2 is larger than that for 
p = due to the increasing repulsive interactions. The 
total energy per particle E/N for an interacting Bose 
gas increases more rapidly than that for an ideal bosonic 
gas due to the finite-size effect, as N increases. There- 
fore, in the Thomas-Fermi limit the kinetic energy of the 
particles in the trap is negligible in comparison to the 
interparticle interaction energy and the trapping poten- 
tial. In addition, Fig. 7 shows the superfluid fraction 
p.s/p along the axis of rotation z as a function of the 
number of particles N at T = 0.4 T* D . The superfluid 
fraction p s /p for an noninteracting gas is also shown in 
the figure. The superfluid fraction p s j p for p, = 0.2 (filled 
circles) is clearly smaller than that for p = (open cir- 
cles). The superfluid fraction p s /p shows no dependence 
on the number of particles N, whereas the interparticle 
interaction plays a significant role in decreasing of the 
superfluid fractions p s /p. Interestingly, when we com- 
pare the superfluid fraction for A = 0.1 to that for A = 
10, we find that the superfluid fraction in the disk-shaped 
trap is significantly smaller than that in the cigar-shaped 
trap because geometric arguments imply that more ex- 
cited longitudinal modes are occupied for A = 0.1 than 
for A = 10. 



B. A disk-shaped trap 

For the disk-shaped trap with the asymmetry factor A 
= u x /u) z = 0.1, i.e. increasing the confinement in the 
axial direction, the length unit is a x = yjfi/muj x and 
the energy unit hw x . In Fig. 5, as we did for the cigar 
geometry, we plot the density profiles n(p) and n(z) as 
a function of p and z, respectively, at T = 0.4 T^ D for 
three different numbers of particles, N — 27 (circles), 64 
(diamonds), and 125 (squares) using p = (open sym- 
bols) and p = 0.2 (filled symbols). It is compared to 
corresponding density profiles n(p) and n(z) for a Bose 
gas with short-range interactions only and for an ideal 



IV. CONCLUSION 

In this paper, Bose-Einstein condensation of trapped 
dipolar Bose gases has been studied using a finite- 
temperature path integral Monte Carlo technique. The 
quantum particles have a dipole-dipole interaction and a 
s-wave scattering interaction. The dipolar potential en- 
ergy is long-range, anisotropic, and partially attractive, 
compared to the short-range and purely repulsive s-wave 
scattering interaction. Using pscudopotential forms to 
describe the interatomic interactions, we have calculated 
the equilibrium properties, such as the energetics and 
structural properties, of a system composed of N dipo- 
lar hard spheres. We find that for a cigar-shaped trap, 
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the net effect of dipolar forces is attractive and shrinking 
of the density profiles is observed; and for a disk-shaped 
trap, the net contribution of dipole-dipole interactions is 
repulsive and the density profiles expand. Accordingly, 
the net contribution of dipolar interactions depends on 
the trapping aspect ratio A. 
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